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I Abstract 

> : 

■ We study BFKL evolution and, in particular, the energy dependence of 

^ . the saturation momentum in the presence of saturation boundaries limiting 

•/^ ! the region of linear BFKL evolution. In the case of fixed coupling evolution 

^ I we confirm the previously found exponential term in Qs{Y) and determine 

O ' the prefactor Y and a dependences. In the running coupling case we find Y^^^ 

r-| ■ corrections to the F^/^ exponential behavior previously known. Geometrical 
scaling of the scattering amplitude is valid in a wide range of momenta for 

Qh! fixed coupling evolution and in a more restricted region for running coupling 

^ I evolution. 

> : 

X 1 Introduction 

In 1983 Gribov, Levin and Ryskin|I|] introduced the idea of parton saturation 
in high energy hard scattering as a dual description of unitarity. Since that 
time our understanding of saturation, and unitarity, in hard reactions has 
progressed considerably 0. We now have a simple model, the McLerran- 
Venugopalan model0, H, ^, which exhibits gluon saturation in a simple and, 
likely, fairly general manner. This model is now being used in order to 
understand general features of heavy ion reactions^, |^, ||, |l^. In deep 
inelastic scattering the Golec-Biernat and WiisthofF model incorporates the 
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essential elements of saturation and gives a surprisingly good fit to much 
of the HERA data for F2 and for diffractive production at small values of 
x]TT|, [T^. There remain, however, many uncertainties in our understanding 
and application of saturation ideas. 

The points we wish to address in this paper are the value and energy 
dependence of the saturation momentum and the form of the scattering am- 
plitude on the perturbative side of the saturation line. To be more specific, 
and in order to illustrate the issues, consider the scattering of a QCD dipole 
of size l/Q on either a hadron or on another dipole of size l//i and with 
relative rapidity Y. Then the saturation momentum, QgiY)^ is the momen- 
tum at which determining the scattering changes from a purely perturbative 
problem, for Q > Qs, to a nonperturbative but weak coupling problem where 
unitarity limits have been reached, for Q < Qs- 

The main question which naturally arises is the F-dependence of Qs and 
what information (dynamics) is actually necessary to control in order to 



calculate the F— dependence of Qs- In general terms one expects BFKL[|13 



dynamics, but not necessarily the BFKL saddle-point solution, to be 



the relevant dynamics since this is the evolution which leads to high density 
partonic systems. Of course one cannot expect linear BFKL evolution to be 
accurate when Q < Qs(Y). As a rough guess one can use the BFKL saddle- 
point solution for high energy scattering at large Q and then define QsiY) 
as the value at which this scattering amplitude reaches its unitarity bound. 
This was done in Ref.l5 for fixed coupling BFKL evolution with the result 
ln(Q2(y)/A2) = ^ where Aq is the solution to x'(Ao)(l - Aq) = 

— x(Ao) with X the usual BFKL eigenvalue function. The danger with this 
procedure is that one cannot justify using the saddle-point method for solving 
asymptotic BFKL evolution even when the scattering amplitude is small. 
(We shall show in Sees. 6 and 7 that this procedure gives the exponential parts 
of the F-dependence of Qs correctly, but misses F-dependent prefactors.) 
This discussion was extended in Ref.l6 to the running coupling case where 
it was also observed that one can expect the scattering amplitude to be a 
function of Q'^/Ql in an extended region where Q'^IQl > 1. In Ref.l7 a 
numerical study of the Kovchegov equation |jl8| was carried out. This has 
the advantage that BFKL dynamics is used when Q/Qs^^ while unitarity 
is imposed in a realistic way when Q/Qs ^ 1- The energy dependence found 
numerically for Qs is close to that expected from simple saddle-point BFKL 
dynamics in the region Q/Qs > 1. 
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In this paper we give a new procedure to solving linear BFKL dynamics 
in the region Q/Q^ > 1 in such a way that the matching with the nonlinear 
dynamics present when Q/Qs ^ 1 should be smooth. To illustrate the main 
problem we face here consider the scattering of a dipole of size 1/Q on a 
dipole of size in fixed (weak) coupling BFKL evolution. Suppose we are 
given QsiY) and that Q/Qs 3> 1. Can we expect the scattering amplitude to 
be given by the saddle-point approximation to BFKL evolution? The answer 
is no, not in general. Because of diffusion there will be many paths, in the 
functional integral sense, which go from /i at F = to Q at F and on the way 
pass through the saturation region. These paths should not be allowed in the 
true solution to the scattering problem. This, however, seems to be the only 
difficulty with using the saddle-point approximation. We face this difficulty 
by converting the usual diffusive behavior in BFKL dynamics to one with an 
absorbing boundary near QgiY), that is we throw away all paths which go 
into the saturation region. 

In order to gain confidence that diffusion with an absorptive barrier is 
the right thing to do, we first study a problem whose answer is known and 
which has many similarities to evolution in the presence of saturation, namely 
non-forward scattering in BFKL evolution. It is well-known that scattering 
at a non-zero momentum transfer q cuts off all infrared dynamics below q. 
In scattering a dipole of size x on a dipole of size xq, with gx, gxo <^ 1 
we expect that one can get the correct, non- forward, BFKL behavior by 
using forward BFKL evolution but with an absorptive boundary 0] which 
eliminates diffusive paths that go into the momenta region below q. In Sees. 
4 and 5 we verify that this is the case. 

In Sec. 6 we evaluate BFKL evolution in the presence of saturation in the 
case of dipole-dipole scattering and where the coupling is fixed. Our main 
results are given in (47) for the scattering amplitude, and in (48) for the 



saturation momentum. The exponent in (48) is as previously found[0, 16 
while the a and F— dependence of the prefactors is new. Eq.(47) exhibits 
geometric scaling|]12| when \n.{Q'^/Ql) > 1, and what is perhaps even more 
remarkable is that there are no unknown prefactors in T. Eq.47 is valid so 
long as ln2(QVQ2) < ^bE^^i [^x^^y, that is within the diffusion regime for 
BFKL evolution in the absence of boundaries. 

In Sec. 7, we deal with the running coupling case. This discussion should 
apply to high-energy hard scattering on protons. Our main results are 
contained in (83), (84) and (85). In (83) we find, for p, = IniQ^A^), 
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the expected leading ^JY behavior |16[ along with a Y^l^ correction while 
(84) and (85) exhibit geometric scaling, but now only for ln(Q^/Q^) < 

[ 7rfc(T-\))x(Ao ) ^] • What is remarkable is that Qg has no knowledge of the 
target, whether that target be a hadron or a small dipole, so long as Y is 
in the asymptotic regime. For a very small dipole target, of size the 

asymptotic regime only begins when Y > 2bN^x{>^o)a'Hf^) ' running cou- 

pling case we are unable to determine target- dependent prefactors in T either 
when the target is a hadron and even when it is a small dipole. 



2 BFKL Evolution; Naive View 

We consider the forward scattering amplitude for a dipole of size 1/Q on a 
dipole of size In the large Nc limit and at the leading logarithmic level, 
this is described by BFKL evolution. When the strong coupling a is fixed, 
the amplitude is 



T{Q,f^,Y) 



dX 



/i2 J 27ri A2(l-A)- 



exp 



2aNr 



TT 



;i-A)ln^ 



where 



X(A)=^(1)--^(A)--^(1-A), (2) 

with ip{X) = r'(A)/r(A) as usual and the integation contour being parallel 
to the imaginary axis with O < Re{X) < 1. F is the relative rapidity of the 
two objects. Expression (1) is symmetric under the interchange Q ^ fi, and 
is normalized at F = to be the elementary cross section of two dipoles at 
the two gluon exchange approximation, which is 



TW.,.y^O)^^(l.ln|). (3) 

In the above, = max((5,/x) and = mm{Q,fj,). 
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Now we define a particular line in the ln((5^/ ii'^)—Y plane by the following 
two conditions 



^"%'(Ao)r + lnf = 0, (4) 



TT li 

and 

^x(Ao)>^-(l-Ao)ln^ = 0. (5) 

Eq.(4) is just the saddle-point condition, while Eq.(5) means that along the 
line = Ql(Y) the exponent in (1) at the saddle-point vanishes. As we shall 
see in a while, this is a line of almost, but not exactly, constant amplitude. 
The solution to (4) and (5) is 

^ = (6) 
x(Ao) 1-A„' ^' 

Ql(Y)=,'e.p(?^^Y). (7) 
The graphical solution to (6) is shown in Fig.l. 



The value of Aq, as determined by (6), is \q = 0.372, which is not too far 
from 1/2, so that F— evolution is the dominant one in the process and the 
approach is reasonable. One can now evaluate the scattering amplitude in a 
region around Ql(Y), by the saddle point method, when aY is large. Since 
Ao is chosen to satisfy the saddle-point condition, we expand x(A) around 
this point. Then by using (4) and (5), the amplitude becomes 



A2(l-Ao)V 



X 



/ 



dX 



2m 



- exp 



l-Ao 



TT 



x"(Ao)y(A-Ao)' + (A-Ao)ln 



(8) 
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Fig.l 



The Gaussian integral is easily done and we obtain 



- 7rin^(gVQg(y)) - 

(9) 

Since higher derivatives of x(A) have been neglected, this result is valid so 
long as I ln(gVg2)| < iaN^x"{>^o)Y/TT. 

A couple of comments need to follow here, (i) When the dipole size 
is inside the diffusion region ln^(g^/QQ) ^ 4Q;x"(-^o)^/7r, and ignoring, for 
the moment, the slowly varying prefactor 1/^/aY, the dominant factor of the 
amplitude is {Qq/Q'^Y~^° ■ This of course has a scaling form, with momentum 
scale Ql{Y) as given by (7); the dependence of this scale is exponential 
and the coefficient in the exponent is known, (ii) The amplitude along the line 

= Ql{Y), as claimed earlier, is close to being constant, but not quite as 
it behaves as 1/ V aY. We are going to find lines of constant amplitude later, 
but an important thing to notice here is, that one can evolve the system to 
arbitrary large rapidity along this line, without facing the problem of making 
the amplitude too big and violating unitarity constraints. This is happening 
because as we increase Y, at the same time we exponentially suppress the 
dipole size. This is in sharp contrast to the usual BFKL evolution, where one 
considers final dipole sizes comparable to the initial one, or more precisely 
inside its diffusion radius ln^((5V/^^) < 4Q;iVcx"(l/2)F/7r. 



T 



A§(l-Ao)V ^/AaN,Y"(Xo)Y 



QliY) 
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3 Diffusion into the Saturation Regime 



The analysis so far is naive, in the sense that we have ignored problems 
arising from the diffusion of the solution, as given in (9), into the satura- 
tion regime. When the forward scattering amplitude becomes of order 27r//x^ 
(or equivalently the amplitude in impact parameter space of order 1), one 
should not trust the solution any more, since unitarity constraints are vio- 
lated. Therefore the dynamics in that region cannot be represented by the 
BFKL evolution. We do not intend to find the amplitude in this saturation 
regime, where presumably we would have to solve the (non-linear) Balitsky- 
Kovchegov equation |]TB|, |T^, but the problem is that the solution in the purely 
perturbative region is altered because of unitarity effects. At this level, one 
can see from (9) that the amplitude becomes of order 27r//i^, when the dipole 
size l/Q becomes \n{Q'^/Ql) ^ — \n{V aY /a^)/(l — Aq). This logarithmic dis- 
tance of from Ql{Y) is small compared to the diffusion radius \/aY when 
aY is large. Thus, even though we start with an initial condition that the 
dipole size is in the perturbative region, as Y increases the diffusion drives 
the dipoles to bigger sizes that enter the saturation regime and therefore it 
will, partly, invalidate the result as given in (9). 

This problem caused by diffusion can be seen from a mathematical point 
of view, when one tries to reproduce (9) by doing two successive evolutions 
in rapidity. Imagine first that we evolve the amplitude from zero rapidity 
to Y/2. Then assume that the amplitude at this rapidity is given by (9), 
with Y ^ Y/2, for dipole sizes 1/Q such that the amplitude is less or equal 
to 27r//i^, and given by 2tt/ix^ for larger dipoles. This serves as an initial 
distribution at y/2 in a simple way to impose unitarity. Then one can evolve 
to find the amplitude at rapidity Y always in the perturbative regime. This 
is a straightforward calculation that we don't present here and one can find 
that the result agrees with (9) provided that V aY < ln(V aY jo?). This is 
of course a small, for our purposes, evolution in Y. It simply states the fact, 
that the solution is incorrect when the rapidity is large enough so that the 
dipoles start diffusing into the saturation regime. 

We shall come back to resolve this issue by imposing unitarity in a more 
proper way in section 6. Before this, and in order to motivate the work 
in that section, we find it useful to consider the non-forward scattering in 
section 4 and how this is related to the diffusion equation in the presence of 
an absorptive barrier in section 5. 
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4 Non-Zero Momentum Transfer 



Just for the purposes of this section and the following one, we adopt a shghtly 

different notation. We consider now the scattering of a dipole of size a; on a 
dipole of size xq, when the momentum transfer of the scattering is q. In this 
case the amplitude is 

(10) 

where 

x(^)=V'(l)-^^(^+^^)-^^(^-^^), (11) 

and 



d^R- 



(|i?+f||i?- 11)1+2^'^ 



(12) 



Following the method presented in [21,22], we calculate the amplitude in the 
Appendix when qx, qxo <S 1. We obtain (for large Y) 



Tnf —Swa'^xxQ 
X < exp 



,{ap-l)Y 



ln^(a;/a;o)^"' 
DY 



— exp 



ln'^{cq'^xxo) 
DY 



(13) 



where ap — 1 = 4aNc{ln2)/7r, D = 56aNcC{3)/7r, the constant appearing is 
c = e^^ /64: and 7 = 0.577.... The result is valid in the domain | ln(a;/a;o)^| <^ 
DY. Here we have done the standard BFKL evolution by expanding the x 
function around its extremum at i/ — (corresponding to A = 1/2 in the 
previous notation). Eq.(13) can also be derived by looking at the scattering 
at a fixed impact parameter b and then taking its Fourier tansformation with 
the integration limits for b being from ~ max(a;,a;o) to ~ 1/q. 
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The first term in the curly bracket of (13) can be recognized as the one 
appearing in the forward case. Notice tliat tlie second term has the same 
diffusion pattern and, because of the minus sign, it will eventually cut off 
large dipole sizes. We will return to comment on this in the next section. 

5 Diffusion in the Presence of an Absorptive 
Barrier 

Going back to (13) and looking at the two terms in the bracket accompanied 
by the 1 / \Jt\DY prefactor, we see that they obey a diffusion equation. What 
is not so clear, for the moment, is the physical interpretation of the second 
term in the diffusion mechanisn, a task that we now turn into. For reasons 
that will be apparent soon we review the diffusion equation in the presence 
of an absorptive barrier, which is 



di 4 df^ 

with the boundary condition 



(14) 



^(-Pa,t)=0. (15) 

We convert this into a Green's function problem, 

/oo 
dp'G(p,p',t-t')V^(p',t'), (16) 
-oo 

where the kernel should satisfy (for p, p' > —pa) 

G(p,p',0)=5(p-p')- (17) 
If we define the Laplace transformation of G{p, p' ,t — t') with respect to time 

by 
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then Eqs.(14) and (17) imply 



(18) 



(19) 



The solution to this equation, satisfying the condition Gi^{—pa,p') — 0, is 
(for p, p' > -pa) 



gUp,p')- 



(20) 



Finally, we perform the integration in (18) along the imaginary axis, to arrive 
at 



G(p, p', t-t')^ 



'■nt 



e *-*' — e *^ 



(21) 



Eq.l6 is quite general and by choosing an initial condition ijj{p',t' = 0) = 
S{p'), which is the relevant one for our purposes, we obtain 



7P{p,t) 



e * — e * 



(22) 



If we let t DY, p \n{xl/x ) and Pa ^ — ln(ca;Qg ), then ip becomes 



ij{x,Y) 



VttDY 





ln^(a;/a;o)^" 




ln^(cg^a;a;o)^ 


|exp 


DY 


— exp 


DY 



, (23) 



which is, apart from the exponential increase and the Sira'^xxo prefactor, 
identical to (13). 

Here we have two equivalent descriptions. Say that we want to calculate 
the forward scattering amplitude with momenta (inverse dipole size) Q 



< 
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q cut-off, where g is a fixed momentum. An intuitive way is to look at 
the non-forward amphtude at momentum transfer q. This is the only new 

scale entering the problem and thus will eventually offer the infrared cut-off 
scale in the problem. An alternative approach is to look directly at forward 
scattering, but at the same time do not allow diffusion, in the dipole size, to 
go into the region Q ^ q. The last must happen in an absorptive way; if the 
dipole hits the boundary at q', it never comes back in the region Q>qaswe 
show in Fig. 2. 



Going back to Eqs.(22), (23) one can expand the exponentials in the 
diffusion region p^,pl^t. Then (22) becomes 

(24) 

Notice the t^^^ power in the denominator. It has the interpretation that the 
probability of not diffusing into Q ^ q,hy cutting out these paths, is 1/aY 
times the probability of all paths. 
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We apply these ideas in the next section, where we return to impose 
unitarity effects in the BFKL evolution considered in Section 2. 



6 BFKL Evolution in the Presence of Satu- 
ration 

Referring to the starting equation (1) for the forward amplitude, we define a 
new hne Ql{Y), which is close to the hne Ql{Y) in the ln((5^///^) —Y plane, 
by the following conditions 



2aNr 



TT 



x'(Ae)r + ln 



01 



(25) 



and 



——x{K)y - (1 - Ac) In — = - ln[ ] . (26) 

TT /J, Z 77 

Eq.(25) is again a saddle-point condition, as Eq.(4), while (26) has an extra 
(3/2) ln(Q;F) term compared to (5). This definition of the critical hne Ql{Y), 
will result in an extra (ccY)^/^ factor in the amplitude, which will cancel 
the (aY)^^"^ factor that is anticipated from the discussion of the previous 
section, more precisely from (24). Therefore we will be able to recover lines 
of constant amplitude. The solution to (25) and (26) is 



(1 - A,)x'(Ae) + x(Ac) = ^4aiy,y (27) 

TT 

and 



Ql{Y) = ^ ^/ ^ - ^ . (28) 
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Eq.(27) determines Ac which is a function of the rapidity Y, and not a pure 
number hke Aq. However, when oY is large, Ac is very close to Aq. In this 
case (27) reduces to 



, , _ 31n[^x--(Ao)y] 

Ac '^°-(i_Ao)^x'W ^ ^ 



In the exponential in (28) we notice that 

x(Ac) _ x(Ao) , 1 x"(Ao) 2, (or.^ 

1^-1^+2 1^^^'^"^'^^ 

as the linear term cancels when we make use of (6). This exponential can 
now be evaluated at Aq, since (29) and (30) imply that the remaining term 
is of order h^{aY)/aY which is small. Then Eq.(28) becomes 

It is clear that the hues Qq{Y) and QliY) are close in the ln((5^///^) — Y 
plane. 

Let's call E the exponent in (1). Expanding x(A) around Ac, this exponent 
takes the form 



E = ^F[x(Ac) + (A - Ac)x'(Ac) + ^(A - Ac)^x"(Ac) + 

TT Z 



-[(1 - Ac) - (A - Ac)] In § - [(1 - Ac) - (A - Ac)] In %. (32) 



Making use of the definitions of Ac and Ql{Y), which are (25) and (26), the 
last expression simplifies to 



^^x"(Ac)nA-Ac)^ + (A-Ac)lng 

2 TT Q't 



13 



3, 



AaN,x"{K)Y' 



(1- Ac)ln 



QV 



(33) 



where the last two terms are independent of the integration variable A. Once 
again we do the Gaussian integration to obtain 



T 



A2(l - Ae) V 



1 3/2 



TT 



01 



l-Ac 



To simplify the notation we define 



exp 



TTln'iQyQl) 



4aN,x"{Xc)Y, 



(34) 



P = In— ,pe = 



and 



(35) 



^ _ 4aiVex"(Ae)y 

TT 

Then the amplitude in terms of p and t becomes 



(36) 



where we defined 



Tra 



A2(l-Ae)V 



-(1-Ac)(p-Pc) 



t^/''lP{p-pc,t), 



(37) 



1 -(P-Pc)^ 

ip{p - pc, t) = —^e t 



(38) 



It is obvious that ip represents the diffusive part of the amplitude since it 
satisfies 
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—il}{p-Pc,t) = -—i){p-pc,t). (39) 

It is at this point that unitarity of the amphtude has to be imposed. 
Eq.(38) is not the proper solution of (39) in the presence of saturation. When 
solving this diffusion equation, we require that we do not include paths which 
go into the saturation region; momenta Q Qs, where Qs is the saturation 
boundary, will be cut out. We do this, following the discussion in the two 
previous sections, by requiring that ijj vanish very close to the saturation 
boundary. Of course in the problem we consider before, the boundary was 
y— independent, since the fixed momentum transfer q was the infrared cutoff. 
Here Qs is F— dependent and it will be a line parallel to the critical one 
Qc- However, Eq.(39) is in terms oi p — Pc and t variables and we can put 
a time independent absorptive barrier in p — p^, which corresponds to a 
F— dependent barrier in momentum Q, thus making our approach for cutting 
momenta Q < Qs{Y) reasonable. 

We still have to exhibit the above in detail and we start by considering 



^s{p, t) = i^ip - Pc, t) - iP{p - pe + 2A, t), (40) 

where is the one given by (38) (the solution with no boundary conditions), 
while in iljs{p,t) we require p > pc — A, as we note that ips = when 
p = Pc — A. A is a parameter that has to be determined. Eq.37 will now 
become 

When p — Pc and A are much smaller than the diffusion radius the 
t-dependent prefactors will cancel with the 1/t factor coming from the ex- 
pansion of the exponentials. In this region wc can also replace Ac by Aq, since 
(Ac — Ao)(p — Pc) <^ (In t)/\/t <^ 1 as implied by (29). Then we are lead to 



T = , 4^(P -Pc + A)e-(^-^o)(^-''^). (42) 

A^(l - Ao) V V^r 
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T 



-2A(a) 




1 



1 -Xo 

Fig.3 



This expression becomes maximum at the point 



p, = Pc-A + -^, (43) 

i — Ao 

which is a finite and a-independent distance away from the point where 
it becomes zero, as shown in Fig.3. Considering the amphtude at p^, we 
determine A by setting 



T{Ps,t)^^, (44) 

with c a constant of order 1 as required by unitarity. Eqs.(42), (43) and (44) 
give 
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We can solve this transcendental equation by iteration, for small coupling a, 
and the solution is 



A(a) = —In — Inln- + 0(const), (46) 

1 — Ao « 1 — Ao a 

where all the constants appearing in the right-hand side of (45) have been 
absorbed in the (irrelevant) constant term in (46). 

We are finally in a position to give a result for the forward scattering 
amplitude, and switching back to our original notation we have 



T 



27r 2a^A{a) 
7v^A2(l-Ao) 



In 



+ A(a) 



(47) 



with Ql{Y) given by (31) and A (a) by (46). 

This result is valid in the diffusion region \n^{Q'^/Ql) <^ Aa N cx" i^o)Y / ''^ 
when Q > Qc (to the right of the critical line in Fig. 4) and in the region 
\\n{Q^/Ql)\ < A(a) when Q < Qc (to the left of the critical line). The 
validity of (47) also requires the rapidity to be large enough, so that 



47r In^a , dir In^iaY) 
V ^ and V ^ ^ ' 

(1 - Ao)2iV,x"(Ao) a 16(1 - Ao)3iV,x"(Ao) a ' 

The first condition is equivalent to 2 A (a) <^ y/t, while the second is a con- 
sequence of replacing Ac by Aq in (28). 

Expression (47) exhibits a scaling behaviour, since the amplitude depends 
only on the ratio Q'^/Ql{Y); lines with constant Q'^/Ql{Y) will be lines of 
constant amplitude. This in fact agrees with recent numerical solutions [p!7 



of the Balitsky-Kovchegov equation, where it was found that the scaling 
behavior is extended in a domain outside the saturation region. 

It is also interesting to notice that, even without the knowledge of the 
exact form of the non-linear effects, one is able to determine the amplitude 
T in (47), with no need to introduce any unknown parameters. 
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Q'(Y) 




Fig.4 



Finally, we come to the issue of the rapidity dependence of the satu- 
ration momentum. In our approach, pg is defined in (43), and by using 
Ps — \n{Ql/ p"^) we are lead to 



Ql{Y)^p^\.jMlJa)a 

where of course an overall constant factor is free. 



(48) 



7 The Running Coupling Case 

Wc now extend our main results, given in (47) and (48), to the case of BFKL 
evolution using a running coupling. Surprisingly this turns out to be not too 
difficult, although we shall not be able to get a result quite as complete as 
(47) in that the overall constant in our amplitude will be undetermined even 
when ln(gVg2) > i. 

Eq.(l) is no longer a good starting point for running coupling BFKL 
evolution. Rather, we write T{Q, p, Y) in the general form 
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T{Q,t,,Y)^a{Q) J^J ^r,^exp[a;y-(l-A)(p-pe)+7 ^Y] 

(49) 

where in this Section p = ln(Q^/A^), not ln((5^//i^), and where pc = Pc{Y) 
and 7 will be specified in a moment. T^^x has no or F— depdendence, but 
does contain the /x-dependence of T. A is the usual QCD A-parameter and 
the a;-integration in (49) goes parallel to the imaginary axis and to the right 
of any singularities T^x may have in u. We now view T as a function of p and 
Y with the //-dependence suppressed. Our normalization is as in the fixed 
coupling case. The BFKL equation is, schematically, 

^lMPl^?^K.iT/a) (50) 

dY TT 

where K is the usual BFKL kernel. Eq.(50) is easily applied to (49) if one 
uses the fact that 

K * e-^^-^^P = 2x{\)e-^^-^^P (51) 

with 

x(A)=V'(l)-^V'(A)-^^^(l-A) (52) 
the usual BFKL eigenvalue function. Using (50) and (51) on (49) gives 

j^T^x exp [uY-{l-X){p- pe) + 7 In y] 

x{. + ^(l-A) + l-?^^^f^,(A)}. (53) 

We now choose pdY) in such a way that when p = Pc{Y) the amplitude T 
becomes almost constant, thus following closely our fixed coupling procedure 
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of Sec. 6. This is done by choosing pdX) so that as much as possible of 
^^^xW is cancelled by (1 — A)^ in the vicinity of the A— values that 
dominate (49) and (53) when p is near Pc(^)- Suppose A— values near Ac are 
dominant. Then write 



X(A) ^ x(Ac) + (A - Xc)x'{K) + ^{X - A.) V(Ac), 



(54) 



in a "diffusion approximation" much like that introduced by Camici and 



Ciafaloni ||23||. Furthermore, when p is not too far from pc write 



a P - 1 

opc Pc 



(55) 



By requiring 



and 



the leading parts of ^^^x(A) are cancelled and the bracket in (53) becomes 



{} 



X 



Tihp 



x(Ae) + (A - A,)x'(Ae) + -(A - A,) V(Ae) 



(58) 



for all terms of the type given in (54) and (55). The final term in (58) is of 
the same form, in (A — Ac) as the second term, but smaller by a factor of 



P-Pc 



so we drop it in what follows. 



If we write 
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T = a{Q) exp[-(l - X,){p - p,) + 7lny]V'(p, Y)To{p^) (59) 

d_ 

dp 



then the factors of (A — Ac) in (58) can be taken to be factors acting on 



ip. Eq.(53) then reads 

which becomes our basic equation. Of course for this whole procedure to work 
the integrals over A in (49) and (53) should be dominated by values where 
A ~ Ac- We can check this at the end by verifying that -^ij) < yi/e'^ so that 
A - Ac ~ Y-^/^ in (49). Before solving (60) we first turn to a determination 
of pc(^) from (56) and (57). 

Comparing (56) and (57) to (25) and (26) one sees that (57) is a saddle- 
point condition while (56) is the condition that T change slowly, depending 
on 7, with Y. Substituting (57) into (56) gives 



x'(Ac)(l-Ac)+x(Ac) = |^ (61) 

analogous to (27). Anticipating that pc{Y)/Y <^ 1 we can expand about Aq, 
satisfying (6), to get 



analogous to (29). 

Now expand (57) about Aq keeping only terms up to first order in Ac — Aq. 
One finds 



Defining po to be the solution to 
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one gets 



with Fq integration constant. Multiplying by pci Eq.(63) can be written as 



_d_ 2 ^ 4iVex(Ao) _ 27P0 

dr^'^ (l-Ao)7r6 (l-Ao)r ^ ' 

where wc have replaced pc by po in evaluating Ac — Aq in (63). Eq.(66) is 
easily integrated to give 



= + const (67) 

1 — Ao 

or 



p,(y) = p,{Y) - (68) 

1 — Ao 

where we have dropped terms of size I/VY on the right-hand side of (68). 
Eq.(68) is the generalization of (31) to the running coupling situation. 

Now we are in a position to solve (60), at least in an approximate way. In 
the denominators in (60) we can replace Pc by po ioi large Y. Then dropping 
Yq in po gives 



^dY Y TY dp^^ ' ^ > 



where 
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' N,{l-Xo){x"{Xo)y 
47r6x(Ao) 



(70) 



and 



c = 



1-Ao 



(71) 



Defining new variables 



and 



Eq.(69) becomes 



C\ 1/3 p - 



yi/6 



t = 6aV3c2/3yi/6^ 



d__d^ 



(72) 



(73) 



(74) 



We do not know how to solve (74) exactly. However, we can find an 
approximate solution of the form 



V' = ^^i(e-A)exp[-^-Ai] 



(75) 



where A is a constant and Ai is the Airy function. One can easily check that 
(75) satisfies (74) up to terms of size ^ip. Since the Airy function decreases 
as exp[— 1^^/^] for large values of ^ our ansatz satisfies (74) in the region 
where i/j is not exponentially small. (Since we wish to use (75) over a wide 
range of values of t it is important that (74) be satisfied including all terms 
of size l/t,^/t and ^^/t as such terms become important after integrating 
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(73) over t. This is indeed the case and this requires the exp[— ^^/t] term in 
(75) whose actual value is small in the region of interest.) 

Return to (59). Using (75) and expressing all variables in terms of ^ and 
t one gets 



have been included with Tq to give Tq. Following our procedure in the fixed 
coupling case we must impose the condition that T not become large in the 
saturation region. Since the first term in (76) grows strongly for negative 
value of ^ we should choose A so that the maximum of T has a value on the 
order of 1///^ (see 44) at ^ = and then goes to zero at ^ = + A, a value 
slightly less than where is the first zero of Ai{^). To that end write 




(76) 



where we have replaced 



l-Ac 
l-Ao 



by 1 in the exponent, and where some constants 



e = A + ei + 5e. 



(77) 



Then 



Ai{^ - A) ~ Ai'{^,)5^ 



(78) 



leading to 



r = exp[-^(ei + 4A)t 




(79) 



T has a maximum at 



Si = 3/t 



(80) 



at which the value of T is 



T 



max 




(81) 
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Clearly we must take 7 = 1 and A = — in order that Tmax be independent 
of F. 

Defining 

p,= (^)'^Vv6[^^ + A + 3/t]+p, (82) 

one finds 

Rewriting (79) in terms of the more physical variables p and Y gives 
or 

when ln(Q7Q2)^_j_ ^ yVe.f'p^ 7-^ and are all related by rather trivially 
calculated constant factors. Eq.(85) is very close to (47), however in the 
present situation we have no control over the value of Tq in contrast to the 
fixed coupling case. 

It was the choice of 7 = 1 in (76) which allowed us to cancel the factor of 
a{Q) ~ present as a prefactor in the right-hand side of (49). That is, 

we have complete control over the Y and Q^-dependence of T near saturation, 
but we have lost our ability to deal with the yU— dependence. Referring to (83) 
it would appear that there are no free factors left in determining Qs, apart 
from an uninteresting possibility of a not large additive constant associated 
with our choice of defining ps at exactly the maximum of T. In particular it 
seems that ps, or Qs, is completely independent of p, in contrast to the fixed 
coupling case where p set the scale Qs- Indeed, we believe this to be the case, 
and this makes the present discussion valid also for the scattering of a dipole 



T' 



(85) 
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of size 1/Q on a proton, however, one must be careful as to the range of Y 
in which (83)- (85) can be apphed. To see the hmits on Y we suppose that 
/i/ A ^ 1 , then for Y not too large there should be a range of Y for which the 
fixed coupling description is applicable even when running coupling effects 
are allowed Referring to (48) one easily sees that so long as 



ln(;.VA^) » ^^i^lMy (86) 

TT 1 - Ao 

Q;(/i) — a{Qs) <^ a(/i), and the fixed coupling approach should be vahd. This 
suggests that the boundary where one must change from a fixed coupling to 
a running coupling description occurs at a transition value of Y 



^ 7r(l - Ao) 1 

If y <^ Ytrans the fixcd coupliug decription is applicable while the running 
coupling description of this section is applicable when Y ^ Ytrans- In order to 
determine Tq in (84), or Tq in (85) we would have to follow the F-evolution 
of the system through the transition rapidity, a task which we have not 
attempted. 

Appendix A 



Here we calculate the amplitude in the non-forward case following pi| , ^ 
With the definitions and notation used in Section 4, we start by exchang- 
ing the denominators of E^^{x) as given in (12), by introducing a Feynman 

parameter (3. After integrating over R we find 

^ d(3 



e-^^-i^2..(V/3(l-/?)gx) (A.l) 

For small xq the exponentials can be approximated by 1 and expanding the 
Bessel function for small argument as 
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K2i. (y(5{l-(5) qx)^^\ r(l + 2iv) 



2iv 



C.C. 



{A.2) 



where c.c. stands for complex conjugate, one finds 



X 



2iu 



c.c.\. 



{A.3) 



The integral over the Feynman parameter (3 is straightforward and we are 
lead to 



r(l + 2iiy) 



r2(l — iu) \qx 



2iu 



qx 



2iu 



(AA) 



The factor outside the bracket in the last expression is just a phase, when u 
is real, and therefore cancels when we consider the product E^'^*{xq)E^^{x), 
for which we have 



X / 



+ V 64 



2iu 



{A.5) 



plus terms odd in u which will vanish when the u integration is done. The 
second term may be written as 



T{l-iiy) 
T(l + w) 



-2^iu 



.27^2 



2iu 



q^xxo \ 
~64 ) 

2'v 2 \ 2ii^ 

e 'q xxq 



3 / V 64 

so that, neglecting terms of order i/^, we finally have 



(A.6) 
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^r(^oX(x) = 2 (^) - 2 {cq'xxof , {A.7) 

with c = Now the u integation in (10) can be done in the standard 

way; by expanding x{^) around 1^ = 0, keeping only up to quadratic terms 
and performing the Gaussian integration. This leads to Eq.(13). The second 
term in (A. 7) gives rise to the second term in (13). 
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